
Fig. 3. — The relative difference in the alignment measures A for the zero and the hrst 
approximations function of S-i . 
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Fig. 1. — The values of T±_ (dotted line) and T|| (solid line) as a function of eccentricity e, 
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The iteration procedure introduced above enables one to calculate (cos 2 f3) and thus 
the measure of alignment with any required accuracy. However, we believe that for majority 
of practically interesting situations no iterations of order higher than one is needed. This is 
not only because a rather high accuracy attainable in the first approximation, but also due 
to more fundamental causes. We adopted the model by Roberge et al. (1993) in which J is 
directed along the axis of major inertia, and this assumption limits our accuracy. In fact, it 
is shown in Lazarian (1994) that the distribution of directions of J in the grain coordinates 
should be described by a Fokker-Planck equation, which solution is sharply peaked only 
when the rotational temperature of a grain is considerably different from that of the grain 
material. Therefore our solution is exact only in the limit of Td/T g — > 0. In our next paper 
(Lazarian, in preparation) we will derive analytics applicable to finite Td/T g . 

Our approach presents several advantages as compared to direct numerical simulations. 
For instance, the formulae obtained above provide a clear physical insight into the 
phenomenon of paramagnetic alignment of interstellar grains; it shows that shapes of grains 
do effect grain alignment, but their influence comes mostly through Si factor. 

By its own, Davis-Greenstein mechanism is not very much interesting for the vast 
regions of the ISM. And it is not only the problem with the magnetic field intensity, 
but mostly with the peculiarities of the polarization curve, which indicate more efficient 
alignment of large grains, while the Davis-Greenstein mechanism tends to go the other way 
round. However, if grains are superferromagnetic, the probability of having a ferromagnetic 
impurity increases with the grain size and the mechanism can be reconciled with the 
polarization curve (see Mathis 1986). The analytical solutions that we obtained, are 
applicable for describing alignment of such superferromagnetic grains. 

This work would not be possible without friendly encouragement by Ethan Vishniac. 
I would like to thank the referee Wayne Roberge for a number of useful suggestions and 
for providing the comparison of my prediction with his numerical results. I am grateful to 
Bruce Draine, Lyman Spitzer and David Williams for illuminating discussions; valuable 
comments by John Mathis and Paul Shapiro are acknowledged. The manuscript was 
improved after my visit to the Princeton University Observatory, and I would like to 
thank Jeremiah Ostriker for providing the financial backing for the visit. The research is 
supported by NASA grant NAG5 2773. 
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heat at constant magnetization and the specific heat at constant field, respectively. The 
difference 



(Caspers 1964), where M is the magnetization of the material, becomes zero for the bias 
field H equal zero and makes x" s i to vanish. In presence of static magnetic field (Morrish 



where r s ; is the spin-lattice relaxation time, // is the Bohr magneton, H s and Hi are, 
respectively, intensities of the static magnetic field and internal field in the paramagnetic 
material; the latter field is ~ Nfi. 

Calculations in Duley (1978) for core-mantle grains provide an enhancement of x" s i 
up to a factor of 10 2 . This is much less that superparamagnetic enhancement could 
provide, but adequate for the purposes of alignment. However, a number of questions 
stay unanswered. Duley studied MgOFeOSiO cores which correspond to the grain models 
in Duley & Millar (1978). Since FeO is antiferromagnetic at least large clusters should 
show no net magnetic moment. Therefore Duley appeals to small particles of size less 
than 100 A, that may in accordance with experiments by Richardson & Milligan (1956), 
Schuele & Deetscreek (1962) and Woods & Fine (1969) demonstrate ferromagnetic ordering. 
Therefore more efficient alignment should be present for the smallest grains, in particular, 
spin-lattice relaxation should occur in bare MgOFeOSiO particles, "since paramagnetic 
defect centers (lattice vacancies, H atoms, other paramagnetic ions) will experience the 
static field from Fe +2 clusters. The resulting enhancement in power dissipation would 
permit 100-200 A MgOFeOSiO cores to be aligned in microgauss fields" (Duley 1978). This 
tendency contradicts to observations which show a better alignment of large grains (see 
Kim & Martin 1994) and therefore the evidence in favor of ferromagnetic relaxation is 
inconclusive. 

It is also important that the ferromagnetic relaxation provides rather limited 
enhancement of the imaginary part of magnetic susceptibility. Therefore grains must stay 
sensitive to magnetic field strength, which makes it possible to use results obtained in 
section 2 to calculate the corresponding measure of grain alignment. 




(63) 



1965) 




(64) 



4. Discussion 
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and find that V p /T should be ~ 1.8 • 10~ 21 cm 3 K -1 , which for T equal 10 K gives 
1.8 • 10~ 20 cm 3 K -1 . Therefore assuming that a unit-sell size for iron is 23.5 • 10~ 24 cm 3 and 
58.9 • 10~ 24 cm 3 for magnetite, they obtain the the optimal number of atoms per precipitate 
particle T which is 1.5 • 10 3 and 7.4 • 10 2 atoms, respectively. Note that according to 
Eq. (58) this is the factor over which the susceptibility Xs P exceeds xo- Therefore n factor 
that enters Eq. (14) increase to < 10~ 7 s -1 for iron and magnetite precipitates as compared 
with 10~ 13 s -1 for ordinary paramagnetic grains; the additional increase in imaginary part 
of susceptibility is due to the fact that this value was chosen to be maximal for frequency 
10 6 s -1 for ordinary materials. This will be true for a very narrow range of T. Indeed, a 
combination 



which controls x'sp changes drastically with V v or equivalently T. In fact, less than 40% 
change in T is enough to reduce n values to those provided by ordinary paramagnetism. 
Provided that the alignment must be rather efficient (<^ 20% corresponds to the data 
in Serkowski et al (1975)), it is unlikely that a substantial fraction of grains contains 
precipitates of exactly required size. Similar criticism of superparamagnetic remedy for the 
alignment problem can be found in Duley (1978). We go further and emphasize that for 
thermally rotating grains, their velocity of rotation scales as /~ 5/ ' 2 ^~ 1/ ' 2 r~ 1/ ' 2 . Therefore we 
do not believe that the Nature fine-tunes to provide grains with precipitates of just required 
size to enable a substantial enhancement of paramagnetic relaxation for a wide range of 
grain shapes, sizes and densities. 

This is even less probable for suprathermally rotating grains. The angular velocity of 
these grains, apart from shape, fractal dimension, and size, depends on the number of sites 
of H 2 formation, and on kinetic energy of H 2 molecules formed over particular types of sites 
(Purcell 1979, Lazarian 1995b). Moreover, tuning is required to be even more fine as for 
rotational velocities close to 10 9 s -1 and the paramagnetic susceptibility of ordinary grains 
reaches its maximal value; therefore it is only T factor of enhancement of n that we are 
speaking about. 

Nevertheless, some enhancement of effective paramagnetic susceptibility is attainable 
for superferromagnetic grains. To understand why this is so, we note, that Davis-Greenstein 
mechanism is based on spin-spin relaxation and ignores the spin-lattice processes (see Jones 
& Spitzer 1967). According to Casimir & Du Pre (1938), the spin-lattice contribution to 
the susceptibility is given by 



where xo(H) is the isothermal susceptibility at a bias field H } Cm and Ch are the magnetic 




(61) 




(62) 



where and cr(°) are the measures of alignment for the first and zero approximations, 
respectively. This function is shown in Fig. 3 for Si in the range from 1 to 10, 7 = 0.2 and 
T m /T a = 0.1. 



3. Enhancement of paramagnetic relaxation 



Superparamagnetic grains were introduced in Jones & Spitzer (1967) in an attempt 
to account for observed alignment through Davis-Greenstein process. Superparamagnetic 
properties of grains emerge, according to Jones & Spitzer (1967), if grains contain 
ferromagnetic cluster within non-magnetic matrix. The susceptibility per volume becomes 

^ 1 T Xo (58) 




UT J 

(Jones & Spitzer 1967), where n is the average number of magnetic precipitates per unit 
volume, fi is the magnetic moment, N is the average number of ion atoms per unit volume, 
T is the average number of iron atoms per a precipitate particle, fi is the "effective" 
magnetic moment of each magnetic ion, Xo i s the ordinary paramagnetic susceptibility for a 
magnetic ion concentration N. 

As the combination Nfi^/Sk is the Curie constant, the susceptibility is enhanced by 
the factor T which is of the order 10 3 — 10 5 . The imaginary part of magnetic susceptibility 
is related to the static value of Xsp(0) m the following way 

where u is the frequency of harmonic oscillations of the magnetizing field, and r is the 
relaxation rate given by 

1 / V P AE\ 

7 = "°»p(-^H- < 60 » 

where the factor v represents the disorienting effects of random thermal modulations which 
non-magnetic matrix produces on a precipitate particle (z/ ~ 5 • 10 9 s -1 for metallic iron 
precipitates), V v is their volume, and AE is the energy barrier restraining the magnetic 
reversal of individual particles (AE/k ~ 6.2 • 10 21 K cm -3 (Brown 1959)). 

If lot < 1, x"p ~ Xs P (0)(ut), while for lot > 1, x % ~ Xs P (0)/(ut). Thus, lot 
determines the relaxation. For a given u, the exponential dependence of r determines a 
rather narrow range of T and this is, to our mind, a serious limitation for the suggested way 
of improving the alignment. For instance, Jones & Spitzer (1967) assume that u = 10 5 s -1 
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Thus 



™' Pi ~ T~l2 I 1 - ^4arcsinyTT^ + 7<Zo^o), (52) 

se n / 



36 

where ^ 

G(ce ) -ae - -ae 2 ,. (53) 

For small ae , G(aeo) tends to 0.257rae o . To determine the behavior of q at this limit, we 
observe that 

xlfa + 1) = 1 + ^-6 1 . (54) 

Therefore for Td/T m Si ^> 1, g ~ (Td/T m Si) 1 <C 1 and the order of the 7-term is smaller 
than in the case Td/T m 6i < 1, ^ > 1, when q tends to unity. In the latter case the order 
of the 7-term is 7ae <C 1. For ae close to unity, the 7-term should be compared with 
the measure of alignment which is of the order of unity, and this ensures a low value of a 
relative error even in the zero approximation. 

Weak alignment 

For ae — > 1, one can use a substitution 



2 _ 1 _ 2 



(55) 



and expand the functions at x = 0. It is easy to see that 

1 / yr^2 _ \ 

cos pi = — 1 arcsmx 

X 2 \ X J 



4 M 

X z \ X 



1 2 1 4 

1 x X 

2 8 . 



X rj X 

x H 1 

6 40 



= ^ + 7jf 2 = ^ + 0.1(l-a$ + 0.17^0, (56) 

where q ~ 0.15(1 — ae 2 ,) = 2SiTd/3T m} which is 0.15(1 — ae 2 ,) for Tj, <C T m . Therefore the 
ratio of the 7-term to the combination (cos 2 p\ — 1/3) that enters the measure of alignment 
is of the order 0.157 <C 1, and this does not exceed 3% even for a maximal value of 7. 

Therefore in both limiting cases of strong and weak alignment the 7-terms are within 
a few percent of the zero approximation. 

To characterize the difference between the zero and the first approximation for arbitrary 
alignment, it is possible to introduce a function 

A = -^5i-- < 57 » 
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It is easy to see that similar equations relate "i" and "i+1" approximations. Indeed, for a 
given cos 2 it is possible to show that 



cos 2 # +1 = — - J 1 - 7g =arcsin^l - ae 2 (l - 7 q 8 ) 1 , (46) 



where ^ 

ft = 2/c , ^ [ cos2 # ~ °- 5 sin2 A ~ + 1) - 1} sin2 A]- ( 47 ) 

Therefore any requested degree of accuracy can be ensured by a sufficient number of 
iterations in accordance with Eqs. (46)-(47). Further we will discuss the limiting cases of 
weak and strong alignment, where some simplifications are possible for the expressions 
obtained. 



2.4. Limiting cases 



The combination 1 — ae 2 , that enters Eqs. (45) and (46) can be rewritten as 



2 r 1 _ T d /T„ 



l-*Z = Si " ■ (48) 
1 + Ol 

Therefore for Tj, <C T m and Si 3> 1, the above expression tends to unity, while for either 
(T m — Td)/T m <C 1 or/and ^ < 1, this expression tends to zero. The first case corresponds 
to strong alignment, the second to weak. Note, that for (T m — Td)/T m <C 1 our expressions 
not directly applicable as the Barnett relaxation is not efficient in such conditions (Lazarian 
1994) and therefore the angular momentum deviates from the direction of the major axis of 
inertia. 

Strong alignment 

If 35q <C 1, we may expand with the accuracy 0{^ 2 ) and 0(^&^) 

1 _ ffi 2 (1 _ 79o) 1 _ Xq ^ ) 

1 - 1 7 , , (50) 



arcsin^l - ae£(l - jq ) arcsin^l - <e 2 Q - 7 ^"u^ - 
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Substituting these coefficients into the Fokker-Planck equation (see Eq. 22) it is easy to 
obtain 

J2 

ln ^ = "2^r m (l- 7 sin 2 /3 ) (35) 

and 

J 2 

ln/j = "2^r a 4l-0.5 7 (l + cos2 / 3 )]' J=x > y (36) 
where, following Jones & Spitzer (1967), we define 

The distribution function 

/ = /,/,/, = const ■ exp [-^^ ^ _ ^ + cos2 A)] ) . (38) 

where 



2 



36^ 



_ T ^h , r m (l + cos^ ) 1 
1 T m \ ' 2T a „(£ 1 + l) J l- 7 sin 2 /V V ^ 

With the adopted accuracy 

ae 2 = ae^l - 79o ), (40) 

where _^ 

9o= an , f1 [cos 2 /3o-0-5sin 2 /3 -{aeg(l + ^)-l}sin 2 /3o]. (41) 

The alignment is independent of J amplitude. Therefore it is possible to perform the 
integration over all possible values of J to obtain 

W x = C\{1 - [1 -ae 2 ]cos 2 /3)- 3/2 , (42) 

where the constant C\ = 0.5aei is defined through normalization 

f sin/3d/3 
Vo (1 - [1 - ae 2 ] cos 2 /3) 3 / 2 1 ' 

Using this angular distribution function it is possible to obtain the first approximation to 
cos 2 (3 

2 W 2 cos 2 f3 sin /3<i/3 

C ° S fh = Jo (l-[l- ffi 2 ]cos 2 /3)3/ 2 - ^ 

which provides 

1 I . seoV 1 _ 7<?o 



cos /3i = — < 1 , arcsim/l — aenfl — 7 <7o) / • (45) 

1 " -o(l " 7<?o) \ ^l- ffi 2 (l- 7 , ) V ° l m \ 1 J 
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where 



2 l + SiT d /T„ 



*° = -"ttjt (27) 



and 

T m = -(T d + T g ). (28) 

The difference of our result as compared to one in Jones & Spitzer (1967) is that we discuss 
alignment of non-spherical grains and therefore our Si is different from the corresponding 
parameter in the latter study; the temperature of grain rotation corresponds to a later 
modification of theory suggested in Purcell & Spitzer (1971) to account for thermolization 
of atoms over the grain surface. In fact, in the same paper we may find interesting 
qualitative considerations relevant to modification of Si for non-spherical grains. Surely, 
Purcell & Spitzer could not get the exact result, as at that time neither alignment of J 
due to the Barnett relaxation, not rapid Larmour precession due to the Barnett induced 
magnetic moment were apprehended, but we can only admire the scientific intuition of the 
two pioneers of grain alignment. The first approximation that we are going to discuss now 
accounts for the non-linearity of the Fokker-Planck equation. 

We will use cos 2 f3 given by Eq. (26) instead of cos 2 f3 in Eqs. (18)-(20). The accuracy 
of the suggested approximation was confirmed by independent calculations in DeGraff et. 
al. (1995). The modified diffusion coefficients can be easily written down: 

f 7 L + tM =7^(1 + ^) (29) 



v gas °mag / °gas 



and 



while 



W(AJ,) 2 > = kT g I b zz (l + l^ 1 + cos2 &)) ' ( 3 °) 

t gas ((AJj) 2 ) mag = 2kT d SiI b zz , j = x,y, (31) 



t gas ((AJ z Y) = kT g I b zz \1 + (1 - 7 sin 2 ft). (32) 



Using Eq. (25), it is possible to write 

WKAJ,-) 2 ) + ((AJ 2 )) ma5 ) = 2kI b zz {T m + T d Si} 

Similarly 



7 (l + cos 2 /3 )r„ 
2(T m + T d Si) 



(33) 



t gas ((AJ z ) 2 ) = kI b z T m (l- 7S m 2 (3 ). (34) 
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where 



<(AJ,) 2 ) = ((AJ.) 2 ), 
((A J,) 2 ) = A(l- 7 sin 2 /3), 



(19) 
(20) 

(21) 



The values ol cos 2 f3 and sin 2 f3 depend on the attained degree ol alignment, and 
therefore to find them, one needs to know the distribution function for grains. In our 
approach, we use the distribution function corresponding to 7 = to obtain the zero 
approximation cos 2 (3 and sin 2 (3 . On the next step, the zero approximation is used to 
obtain the first approximation for the distribution function. This process may be continued 
until the required accuracy is attained. 



2.3. Iteration procedure 



We start with assuming 7 = 0. In this case, the separation of variables is applicable to 
Eq. (3). The stationary equation for the z component is as follows 

\w^ {AJz)2)fz) ~ w z i{AJz)fz) = 0> (22) 

which gives a formal solution 

ta/ * = -wfe> + comt '- (23) 

In the case of x and y components, one needs to account for paramagnetic relaxation and 
therefore 

J 3 ( 1 1 \ 

((AJj) 2 ) + ((AJj) 2 ) m(l5 \~tgas tmag J 

To characterize the relative importance of magnetic torque, it is advantageous to define 

^ _ (A Jg) mag _ tgas _ 3 K,V B 2 

1 (AJj) t mag 40? nmvthbUw ' 

With thus defined 8\ the problem is similar to that discussed in Jones & Spitzer (1967). We 
refer our reader to this work and write here the answer 



2 /9 1 

cos po 



1 — aeg 



1 , ° arcsin\/l — ae, 2 

1 — ae 



(26) 
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where the subscript "i" stands for x,y and z while 

O Jb 

t = — — (11) 

The diffusion coefficients discussed above characterize gas-grain interactions. The effect 
of magnetic field on grains is imprinted through the coefficients 

(AJ 3 ) mag = j = x,y (12) 

mag 

(AJ z ) mag = (13) 

(see Jones & Spitzer 1967) and 



I 

kVB 

where for slow rotation, k 2.5 • fO -12 ^ -1 s (Spitzer 1978). Similarly 



t = — ^— (14) 

b mag T/ D2 ' V / 



{(AJj) 2 )^ = 2kT d VB\, j=x,y (15) 

and 

((AJ,) 2 ) ma , = 0, (16) 

where Eqs (13) and (16) reflect the fact that magnetic field does not influence grains 
rotating about its direction. 



2.2. Small parameter 



It is a common knowledge that asymptotic analytical results are usually attainable if 
a small parameter is present in a model. The difficulty in direct solving the Fokker-Planck 
equation stems from the fact that ((A J;) 2 ) depends on the direction of J in respect to 
the magnetic field. Jones & Spitzer (1967) did not face this problem as this dependence 
disappeared for spherical grains for which T|| = 

In general, T|| differs from but according to Fig. 2 their ratio is close to unity and 
therefore the parameter 7 

7 = 1 ~ (17) 

does not exceed 0.2. We choose it to be the small parameter of our model. Thus Eqs. (4)-(6) 
may be rewritten as 

((AJ,) 2 ) = A(l - 0.5 7 [1 + cos 2 /3]), (18) 
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where / is a distribution function of the angular momentum J, while (A J;) and (A J; A Jj) 
are coefficients. An important job of calculating the above coefficients was done in Roberge 
et al. (1993). Wherever possible we preserve the notations adopted in the said study to 
facilitate a comparison of our predictions with the numerical results obtained in Roberge et 
al. (1993). In fact, such a comparison was done in DeGraff, Roberge & Flaherty (1995) and 
proved an exellent agreement of the results obtained through these two approaches. 



It is shown in Roberge et al. (1993) that only the tensor components 
<(AJ,) 2 ) = ^nmbivl (l + ^ [I\(l + cos 2 fi) + r„ sin 2 



(4) 



((AJ y Y) = ((AJ X ) 2 ) (5) 

and 

((A J,) 2 ) = ^nm^blvl (l + |) (I\ sin 2 /3 + I,, cos 2 /?) (6) 
are not equal to zero after accounting for the Larmour precession. 2 Above Tj, and T g are, 



respectively, dust and gas temperatures, v t h = y2kT g /m is the thermal velocity of gaseous 
atoms of concentration n and mass ra; the coefficients Tj_(e m ) and rj|(e m ) are geometrical 
factors 

Tx(e m ) = ^{7 - e 2 m + (1 - e 2 Jg m (e m ) + (1 - 2e 2 J[l + e" 2 (l - [1 - e m ] 2 5 (e m ))]}, (7) 
ri|(e m ) = ^{3 + 4(1 - e 2 Jg m (e m ) - e 2 m [l - (1 - e^) 2 g(e m )]} , (8) 



16 



wi 



ith 



g(e m ) = —ln( 1 -^). (9) 



2e m \ 1 



(see Fig 1) 



Note, that r_i_(e m ) and F||(e m ), coincide for both spherical grains (e m — > 0) when they 
attain the value of 1, and flakes (e m — > 1), when they attain the value of 3/8. The fact that 
cos f3 enters the expression for the diffusion coefficients in Eqs (4) and (6) reflects the fact 
that grain disorientation due to gaseous bombardment depends on grain geometry. 

As for the mean angular momentum increments, it was obtained in Roberge et al. 
(1993) 

(AJ 8 ) = (10) 

£ gas 



2 This precession takes place at a scale ~ 10 5 s if the Barnett effect endows a rotating grain 
with a magnetic moment. This time is substantially smaller than that of grain alignment. 
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to the Barnett relaxation, vector J is directed along the major axis of inertia 1 , which for 
an oblate spheroid coincides with its short axis. Therefore within the adopted model the 
measure of J alignment coincides for oblate grains with the Rayleigh reduction factor. 

Davis-Greenstein mechanism aligns grains rotating at thermal velocities. Such a 
rotation is expected for grains in molecular clouds, where the concentration of atomic 
hydrogen is low. Indeed, according to Purcell (1979), recoils from nascent H 2 molecules are 
likely to be the major driving force for suprathermal rotation. It is claimed in Lazarian 
(1995a) that aromatic carbonaceous & graphite grains do not attain suprathermal angular 
velocities due to H 2 formation. Therefore Davis-Greenstein process may also be applicable 
to some fraction of grains within the diffuse ISM. 



2. Davis-Greenstein process 



2.1. Diffusion coefficients for a model grain 



For our present purposes, we will consider only oblate spheroidal grains as there are 
indications that aligned grains are oblate rather than prolate (Aitken et al. 1985, Lee & 
Draine 1985, Hildebrand 1988). The spheroid semiaxes parallel and perpendicular to the 
grain symmetry axis are denoted a m and b m ib m > a m ), respectively. The core-mantle 
interface is assumed to be a spheroid confocal with the mantle surface. The core semiaxes 
are denoted by a c and b c . The eccentricity of the core (i = c) as well as that of the mantle 
(i = m) is 



and may be different for different components of the grain. 

It was shown in Jones & Spitzer (1967) that paramagnetic alignment can be described 
through the Fokker-Planck equation (see Reichl 1980) 



1 A quantitative discussion of Barnett relaxation in Lazarian(1994) showed that this 
approximation is true with high accuracy only when the grain rotational temperature 
substantially exceeds that corresponding to grain material. 
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1. Introduction 

Discovered nearly half a century ago (Hiltner 1949, Hall 1949) the alignment of the 
ISM grains remains a long standing unsolved problem (see Goodman et al 1995). To solve 
this problem is not only challenging, but alluring, as this should enable one to use the 
wealth of polarimetric data for quantitative studies of the ISM. 

Paramagnetic alignment of thermally rotating grains is known as the Davis-Greenstein 
process. Introduced as far back as 1951, this process was later criticized for not being strong 
enough (Spitzer 1978, Whittet 1992). To improve the mechanism its superparamagnetic 
and ferromagnetic modifications were introduced (Jones & Spitzer 1967, Duley 1978, 
Mathis 1986). At present, however, observational evidence for paramagnetic alignment is 
inconclusive (Hildebrand 1988) and we believe that one of the reasons for this unsatisfactory 
state of affairs has been the absence of the analytical theory of alignment for non-spherical 
grains. Indeed, alignment of non-spherical grains was studied through Monte-Carlo 
simulations (see Purcell 1969, Purcell & Spitzer 1971), through numerical solving of 
Langevin equation in Roberge, DeGraff & Flaherty (1993), but analytical results were 
obtained only for spherical grains (Jones & Spitzer 1967). 

In present paper we address this problem and propose an analytical approach 
to studying alignment of non-spherical grains. First (section 2), we formulate the 
Fokker-Planck equation for non-spherical grains and find a small parameter inherent to 
the problem. Then we use a perturbative approach to solve the problem analytically. The 
recursive formula that we obtain enables one to calculate the measure of alignment with 
any requested degree of accuracy. However, we show, that even the first approximation is 
accurate within 1%, this value that should suffice present-day requirements. Further on 
(section 3) we show, that the assumption of superparamagnetism requires fine tuning of 
sizes of superparamagnetic particles, and therefore is improbable, but the assumption of 
superferromagnetic inclusions is more difficult to dismiss. The ferromagnetic modification 
of the theory corresponds to the two orders of magnitude increase in values of imaginary 
part of paramagnetic susceptibility and, as a result, superferromagnetic grains should be 
aligned in accordance with the analytical theory introduced in this paper. 

If grain axis makes angle f3 with magnetic field H, for an ensemble of spheroidal grains 
the measure of alignment can be described by the Rayleigh reduction factor (Greenberg, 
1968): 

^<cos 2 /3-i), (1) 

where here and further on angular brackets () denote an ensemble averaging. In the model 
of Roberge et al (1993) that we also adopt for the rest of the paper, it is assumed, that due 
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ABSTRACT 

Paramagnetic alignment of non-spherical dust grains rotating at thermal 
velocities is studied. The analytical solution is found for the alignment measure 
of oblate grains. Perturbative approach is used for solving the problem. It is 
shown that even the first approximation of the suggested iteration procedure 
provides the accuracy well within one percent of the expected measure of 
alignment. The results obtained are applicable both to paramagnetic and to 
superferromagnetic grains. 
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